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Abstract. In application to positron annihilation spectroscopy, Ps atom is considered not as a 
point particle, but as a finite size e + e~ pair localized in a bubble-state in a medium. Variation 
of the internal Coulombic e + -e~ attraction vs. the bubble radius is estimated. 

Introduction Typical lifetimes (up to annihilation) of a para-positronium atom(p-Ps; spin = 
OJj in condensed medium are about 130-180 ps. They are close to the p-Ps lifetime in vacuum 
(125 ps). The ortho-positronium lifetime in a medium is considerably shorter (about 100 times; 
some ns) in comparison with that in vacuum. This is due to the so-called pick-off process - 
prompt 27-annihilation of the e + , composing Ps atom, with one of the nearest e~ of surrounding 
molecules, whose spin is antiparallel to the e + spin. Just this property turns Ps into a nanoscale 
structural probe of matter. The theoretical task consists in calculating the pick-off annihilation 
rate X po , i.e. in relating X po with such properties of the medium like surface tension, viscosity, 
external pressure and size of the Ps trap. 

Originally, to explain the unexpectedly long lifetime of the ortho-Ps atom in liquid helium 

i/-) . R.Ferrel pQ suggested that the Ps atom forms a nanobubble around itself. This is caused by a 
strong exchange repulsion between the o-Ps electron and electrons of the surrounding He atoms. 

i/-) I Ferrel approximated this repulsion by a spherically symmetric potential barrier of radius R^. 

^ I To estimate the equilibrium radius of the Ps bubble he minimized the sum of the Ps energy in 
a spherically symmetric potential well, i.e. ir 2 h 2 /AmR 2 ^ = -^(vras/i^oo) 2 , Ry=13.6 eV, and the 
surface energy, 47ri? 2 XD o", where a is the macroscopic surface tension coefficient. The following 
relationship is hereby obtained for the equilibrium radius of the bubble: 
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The Tao-Eldrup model Ferrel's idea got further development in the studies of Tao [2] 
and Eldrup et al. [3]. They considered the Ps atom as a point particle in a liquid, i.e. in a 
structureless continuum, Fig. ??. The repulsive Ps-liquid interaction was approximated by a 
rectangular infinitely deep spherically symmetric potential well of radius R^. In such a well, 
the wave function of a point particle has the following standard expression: 

*(0 < r < R^) = sin (WjU ^ ^ (r > Roo) = Q (2) 

^/2nR 00 r 

Here, r is the Ps center-of-mass coordinate. Because the Ps wave function equals to zero at 
the bubble radius (and outside), there is no e + overlapping with outer electrons of a medium. 
So, pick-off annihilation is absent. To overcome this difficulty it was postulated that molecular 



1 Ps is a bound state of e + and c 
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Fig. 1: Dependence of the pick-off annihila- 
tion rates jl] vs. surface tension in differ- 
ent liquids. Solid curve shows the correlation 
given by the Tao-Eldrup at A + = 2 ns _1 and 
optimal value 6 = 1.66 A (obtained from 
fitting of these data by means of Eq. (J3])). 
Dashed curve illustrates simplest approxima- 



tion \ po OC rr 



1/2 



electrons, which form a "wall" of the Ps bubble, may penetrate inside the potential well. This 
results in the appearance of a surface layer of thickness 5 = -R^ — R having the same average 
electron density as in the bulk. As a result, the pick-off annihilation rate X po becomes non-zero. 
It is proportional to the e + overlapping integral with the electrons inside the bubble: 
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This is the well-known Tao-Eldrup formula. Here, A + ~ 2 ns _1 is the e + annihilation rate in an 
unperturbed medium (it is proportional to Dirac's 27-annihilation cross-section and the number 
density of valence electrons). The thickness 5 of the electron layer is an empirical parameter, 
which may have different values in various media. Substituting Eq. (TTJ) for R™ into Eq. (J3J), one 
obtains the relationship between X po and a with one adjustable parameter, S. g It may be easily 
obtained by fitting experimental pick-off annihilation rates with the relationship ([3]), Fig. [TJ 
Thus we obtain 5 ~ 1.66 A. Eq. ([3]) with this value of 5 is widely used for recalculation of the 
observed pick-off annihilation rate into the free volume 47ri? 3 /3 of the cavity, where Ps atom 
resides and annihilates. 



Further development of the Ps bubble models Along with the development of the 
"infinite potential well" Ps bubble model, another approach based on the finite potential well 
approximation was also elaborated [SJ El El HJ- However in both approaches, the Ps atom 
was approximated by a point particle. This leads to a significant simplification, but it is not 
justified from a physical viewpoint, because: 

1) the size of the localized state of Ps (size of the Ps bubble) does not significantly exceed 
the distance between e + and e~ in Ps; 



2 If we make use of the relationship 
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Thus we obtain X po ex c 1 / 2 , which is displayed in Fig. [TJ 
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2) during the formation of the Ps bubble there is a substantial variation of the Ps internal 
energy (particularly of the Coulombic attraction of e + and e~), which is completely ignored in 
the "point-like" Ps models. In a vacuum or in a large bubble, the internal energy of Ps tends 
to — Ry/2 = —6.8 eV. In a continuous liquid (no bubble) with the high-frequency dielectric 
permittivity e ~ n 2 (n rs 2-3 is the refractive index) the energy of the Coulombic attraction 
between e + and e~ decreases in absolute value by a factor e 2 ra 4-9. The same takes place 
with the total Ps binding energy, which tends to the value — Ry/2e 2 rs —(1-1.7) eV (this is a 
simple consequence of the scaling e 2 — >■ e 2 /e of the Schrodinger equation for Ps atom). Thus, 
the change in the Ps internal energy during Ps formation may reach 5 eV. Obviously, this 
represents an important contribution to the energetics of Ps formation. The aim of the present 
work is towards a more accurate estimation of this contribution, that has not been done yet. 

There is only a small number of papers where the consequences of the finite size of Ps are 
discussed in application to positron annihilation spectroscopy. To calculate A po , the Kolkata 
group [H] suggested to smear the Ps atom over the relative e + -e~ coordinate exactly in the 
same way as it is in a vacuum. Such an approach is valid for rather large bubbles. However, 
they do not discuss the variation of the internal Ps energy. 

In [10] the Ps atom is considered as a finite sized e~e + pair, but the variation of the 
Coulombic interaction because of dielectric screening is not discussed. It was assumed that 
e~ is confined in an infinite potential well and e + is bound to it by means of the Coulombic 
attraction. The wave function of the pair was taken as a series of orthogonal polynomials, their 
weights being determined from a minimization procedure of the total energy of the pair. 

Hamiltonian of e + e pair in a medium Let the e + e~ pair (Ps atom) have already formed 
in a liquid a nanobubble (spherical cavity; Ps bubble) of radius R (the onset of coordinates is 
taken at the center of the bubble, Fig. [2]). Together with the molecules surrounding the e + e~ 
pair, one has to deal with a quite intricate many-body problem with a complex hamiltonian. 
We reduce it to the following form: 

H « - ^ (A 2 + ^ A ' } + U(r + ) + U(r_.) - U c (t + , r_, R, e). (4) 

Terms with Laplacians A + and A_ over r + and r_ (e + and e~ coordi- 
nates) stand for the kinetic energies of the particles. U(r + ) and C/(r_) 
describe the individual interaction of e + and e~ with the medium. For 
them we adopt the following approximation: 
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Here, V + and V ~ are the e + and e~ work functions, respectively (Vq is 

a commoner notation for the electron work function). The work function is usually introduced 

as the energy needed for an excess particle to enter the liquid without any rearrangement of 

its molecules and to stay there in a delocalized state, having no preferential location in a bulk. 

One may say that V + and V ~ are the ground state energies of the quasifree e + and e~ , because 

their energies at rest after having been removed from the liquid to infinity are defined to be 

zero. 

Vq consists of 1) the e~ kinetic energy, arising from its exchange repulsion from the "core" 

electrons of molecules (atoms), and 2) the energy due to the polarization interaction of e~ with 

the mediumo According to the theory of the quasifree electron [11] , this polarization interaction 




3 In case of e + the kinetic contribution to Vq~ is due to the Coulombic repulsion from the nuclei (the exchange 
repulsion is absent). 



Table 1: Electron work function for different liquids at room temperature [12 



Liquid 


V ~, eV 


Liquid 


V ~, eV 


helium; 4.2 K 


1.3 


benzene 


-0.14 


n-dodecane 


0.2 


isooctane 


-0.17 


n-decane 


0.18 


toluene 


-0.22 


n-heptane 


0.12 


neopentane 


-0.38 


n-hexane 


0.1 


MeOH, EtOH, PrOH 


-0.4 


nitrogen; 77.3 K 


0.05 


xenon; 170 K 


-0.57 


n-pentane, c-hexane 


0.01 


water 


-1.2 


argon; 86.4 K 










may be estimated as a sum of two parts: a) interaction of the e~ with the molecule where it 
resides, £/l nt (to calculate [7i nt the electron is considered as an electron cloud smeared over the 
molecule), and b) interaction of the e~ with all the other molecules, C/f* = (1 — l/e)e 2 /2i?ws, 
(this expression is similar to the well-known Born formula for the electron solvation energy). 

Experimental values for V Q ~ are known for many liquids (Table [T]). Because of a lack of 
experimental data on the e + work functions, we shall admit that they are approximately the 
same as for e~: V Q + ~ V Q ~. So we may conclude that \V + + V ~\ < 1 eV. Thus \V Q + + V ~~\ is less 
than the variation of the internal energy of the pair, ps Ry(l — l/e 2 )/2 rj 5 eV, related with the 
variation in the dielectric screening of the e + -e~ attraction in the bubble formation process. 

Note that the use of Eqs. (jSJ) for the potential energies of the e + and e~ interaction with 
the medium, assumes that the polarization interaction remains the same whether e + and e~ 
(both in the quasi-free states) are well separated or form the quasi-free Ps atom. Since, for 
distances larger than the size of a molecule, qf-Ps is nearly an electrically neutral particle, the 
contributions UZ ut ~ U™*, which come from a long-range polarization interaction of the quasi- 
free e + and e~ with the medium, should be absent in U(rJ) + U(r + ) in Eq. (J41). Therefore, it is 
reasonable to consider at least two cases: 1) when the above mentioned polarization correction 
is neglected and V ~ + V + — > and 2) when the terms U° ut « U™ 1 ~ — 1 eV are subtracted 
from the work functions and V ~ + V + — > 2 eV. Both cases are considered below. 

In Eq. (jlj) U c stands for the Coulombic interaction between e + and e~ in a polarizable 
medium. Assuming that the medium has the dielectric permittivity e of the bulk and a spherical 
cavity of radius R (inside the cavity e — 1), one may calculate U c by solving the Poisson 
equation. Denoting the e + and e~ coordinates as r + and r_, U c may be written in the form of 
the following series via the Legendre polynomials Pi(x = cos 9) [13J: 
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Fig. 3: Dependence of the e + -e~ Coulombic interaction energy for different locations of e + and 
e~ around the bubble (here we adopt the radius of the bubble R equal to Aa B w2 A). z + and 
2_ are the e + and e~ displacements from the center of the bubble along the z-axis. The dashed 
curves describe the unscreened (red) and completely screened (e times less; blue) Coulombic 
energies between e + and e~. 



Here, the argument of the Legendre polynomials is x = cos??, where •& is the angle between the 
z axis and the direction of r„. Note that the summation of these series is simplified considerably 
when using the following recurrent relationship 

P,(z) = [(2/ - 1)xP w (i) - (/ - l)P t . 2 (x)}/l. 

Particular dependencies of U c for some selected arrangements of e + and e~ and the cavity are 
shown in Fig. [3j Thus, we are able to take into account the variation of the e + e~ Coulombic 
energy during the formation of the Ps bubble. Similarly, the dielectric screening is used in the 
polaron problem and the ion-electron recombination problem (Onsager's formula) [Til H3] . 



Wave function of the e + e _ pair and minimization of its total energy (H) Keeping in 
mind further use of the variational procedure, let us choose the normalized e + e~ wave function 
in the following simplest form: 
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In both cases of a rather large bubble and a uniform dielectric continuum, \l/+- breaks into a 
product of two terms: the first one depends on the distance r between e + and e~, and the second 
one depends on the center-of-mass coordinate r cm . Parameters a and b are the variational ones, 
over which we have minimized the energy of the e + e~ pair: 



E(a,b,R) = ($ + -\H\ty + _) ->■ min 



a(R), b(R). 



(8) 



The simplest verification of the calculations is to recover two limiting cases. In case of large 
bubbles (R — > oo), one should reproduce the "vacuum" state of the Ps atom: its total energy 
must tend to — Ry/2 = —6.8 eV, the kinetic energy to +Ry/2 and the Coulombic energy 
to — Ry. In case of small bubbles (R — > 0), the delocalized qf-Ps state must be reproduced. 
The Schrodinger equation for qf-Ps has the same form as for the vacuum Ps, but with the 
substitution e 2 — > e 2 /e. Then the total qf-Ps energy tends to V$ + V^ - — Ry/2e 2 , its kinetic 
part tends to +Ry/2£ 2 = 1.7 eV (e = 2) and the Coulombic energy tends to — Ry/e 2 = —3.4 
eV. Fig. H] displays optimal values of a and b as well as different contributions to the total 
energy of the e + e~ pair when Vq + Vq = and 2 eV. 
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Fig. 4: Dependencies of the optimal parameters a and b vs. R, the bubble radius. They enter 
the e + e~ wave function and yield the minimum of the total energy (H). The relative contact 
density r] c and different energy contributions to (H) (at optimal a and b) are shown as well. 
The upper drawings correspond to the case V + + V- = and the lower ones to V + + V- = 2 
eV. In both cases it was assumed that e — 2. 



Relative contact density and pick-off annihilation rate In the framework of the devel- 
oped scheme using the wave function ([7]) it is easy to obtain the relative contact density r\ c in 

Ps atom: 

f ftZ 3 r + rf 3 r_|fr + _(r + ,r_)| 2 (5(r + -r_) g\ 

Vc J/rf 3 r + d 3 r_|^^(r + ,r_)| 2 5(r + -r_) " a 3 (R)' U 

This quantity determines the observable Ps annihilation rate constant (including the case with 
applied permanent magnetic field). The resulting dependencies of r] c are shown in Fig. H] (on 
the left). Because, for qf-Ps, parameter a is equal to eclb, for qf-Ps the value of r\ c should be 
1/e 3 = 1/8, which is well recovered in numerical calculations. When R increases, r] c approaches 
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Fig. 5: Pick-off annihilation rate constant of Ps, lo- 
calized in a bubble of R when V+ + V- = (green 
curve) and V + + V- = 2 eV (brown curve). For 
small R (< 2 A) the calculated values of X po are 
equal to A + . The red line shows pick-off annihi- 
lation rate constant, calculated according to the 
Tao-Eldrup formula. The dashed line is the calcu- 
lation according to the finite potential well model 
(for comparison we adopted that the depth of the 
well is (1 — l/e 2 )Ry/2 ^5.1 eV and its radius is R. 
The minimal radius of the well when there appears 
an energy level is 5jj = 0.96 A). 



unity, because a tends to its vacuum value as- Knowing the expression for the wave function 
([7]), one may calculate the positron overlapping Pr with molecular electrons, surrounding the 
Ps atom, and therefore find out the pick-off annihilation rate constant: 



Xpo(R) ~ A+Pr, 



r + >R 



d r_ |\Eu_(r_L, r_ 
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Here, A + ~ 2 ns _1 is the annihilation rate constant of "free" positrons. Results of calculations 
of X po (R) for optimal a and b values, which correspond to the minimal Ps energy at a given R, 
are shown in Fig. [51 



Results and discussion 1) It is usually considered that Ps is a solvophobic particle, i.e., 
it forms a bubble when entering a liquid because of exchange repulsion between e~ in Ps 
and the surrounding molecular electrons. If the work functions of e + and e~ are negative 
(Vj} + ~ Vq' < 0), each particle considers a cavity as a potential barrier. So they are pulled to 
the bulk by polarization interaction with the medium. Nevertheless, even in this case the Ps 
bubble may be formed due to an enhancement of the Coulombic e + e~ attraction inside the 
cavity (no dielectric screening inside). This feature cannot be taken into account when Ps is 
simulated as a point particle. 

2) It is seen that the behavior of the total energy of the pair (red curves in Fig. HJ) strongly 
differs from the Tao-Eldrup prediction (green dashed curves; the first term in Eq. (JTJ), where 
Roo is replaced by R), as well as from the expectation based on the finite potential well model 
(brown curves in Fig. HJ the Coulombic potential cannot be approximated well by a rectangular 
spherically symmetric potential). The same is true for the pick-off annihilation rate, Fig. 

3) Calculations demonstrate one common feature: up to R < 1.5 — 2.2 A all dependencies 
remain the same as in a medium without any cavity, but at larger R there are significant 
deviations. This is related to the known quantum mechanical phenomenon - absence of a bound 
state of a particle in a small finite 3d-potential well. In such cavities, Ps cannot be bound, it 
does not exert any repulsive pressure on their walls and does not stimulate their transformation 
towards the equilibrium Ps bubble. The possibility of finding a suitable preexisting cavity, 
sufficient at least for preliminary localization of qf-Ps, may be a limiting factor for the formation 
of the Ps bubble state. 

4) One may find an equilibrium Ps bubble radius by minimizing the sum of the total e + e~ 
energy (H) and the surface energy of the bubble. For water it turns out to be 5-5.2 A which 



is about 2 A larger than predicted by the Tao-Eldrup model. For such a large bubble, the 
relative contact density is r] c « 0.9, Fig. HI It is somewhat higher than the experimental values 
(0.65-0.75 [IS])- This discrepancy may indicate that e + and e~ really interact with a medium 
in a different way, for example, Vq > V^ . It means that the Ps electron may be trapped by 
a cavity, and e + will be bound to this trapped e~ by the Coulombic attraction. This scenario 
may be also considered in the framework of the present approach, but the expression for the 
trial wave function of the new pair must be written in an "asymmetric" (towards e + and e~) 
form: 

* + _(r + , r_) « eXpHr+ - r -^- r - /2&) . (11) 

87rva 3 6 3 

5) Any Ps bubble model reduces the original many-body (multi-particle) problem to a 

simpler one, that of one or two particles in an external field, which simulates the interaction with 

the medium. To calculate this field one usually relies on some macroscopic approaches. However, 

their validity always remains uncertain (for example, how to relate the actual arrangement of 

molecules around the Ps bubble with the jump of dielectric permittivity outside the bubble and 

so on). 
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